Dephasing of solid-state qubits at optimal points 
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Motivated by recent experiments with Josephson-j unction circuits, we analyze the influence of 
various noise sources on the dynamics of two-level systems at optimal operation points where the 
linear coupling to low-frequency fluctuations is suppressed. We study the decoherence due to non- 
linear (quadratic) coupling, focusing on the experimentally relevant 1/f and Ohmic noise power 
spectra. For 1/f noise strong higher-order effects influence the evolution. 
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For quantum-information processing it is crucial to 
reserve phase coherence. Recent experiments III 00 [J 
with Josephson-junction circuits demonstrated long- 
lived coherent oscillations. They showed resolution, suf- 
ficient for detailed studies of the dephasing times and de- 
cay laws, stressing the need for the theory analysis of the 
dissipative dynamics of qubits subject to relevant noise 
sources. 

In solid-state systems decoherence is potentially strong 
due the host of microscopic modes. In Josephson qubits 
the noise is dominated by material-dependent sources, 
such as background-charge fluctuations or variations of 
critical currents and magnetic fields, with power spec- 
trum peaked at low frequencies, often 1/f. A further 
relevant contribution is the electromagnetic noise of the 
control circuit, typically Ohmic at low frequencies. The 
1/f noise appears difficult to suppress and, since the de- 
phasing is dominated by low-frequency noise, it is par- 
ticularly destructive. On the other hand, Vion et al. |2j 
showed that the effect of this noise can be substantially 
reduced by tuning the linear longitudinal qubit-noise cou- 
pling to zero. The same strategy, suppressing the linear 
qubit-detector coupling, was used to minimize the effect 
of the quantum detector in the off-state. The achieved 
coherence time was 2-3 orders of magnitude longer than 
in earlier experiments. 

The 1/f noise (more generally, strong low-frequency 
noise) plays a major role in many solid-state systems. 
The long-range correlations in time make the analysis of 
its effect difficult. In this letter we analyze the dynamics 
of a qubit subject to singular noise, with a focus on non- 
linear coupling and higher-order effects. We also apply 
the developed formalism to an environment with regular 
(e.g., Ohmic thermal) noise. 

The dynamics of a dissipative two-level system (spin 
1/2, qubit) can be described by the Hamiltonian: 

H = -~ (e a z +V a z + U a x ) + H hath ■ (1) 

The longitudinal (V) and transverse (U) fluctuations 
may result from various microscopic noise sources, de- 
scribed by flbath- I n this letter we analyze the dephasing 



in the situation of quadratic longitudinal coupling, 

V = XX 2 , U = 0, (2) 

to a source of Gaussian noise X(t), with noise power 
Sx{w). Here X{t) is a basic physical quantity (e.g., 
voltage or magnetic field), which controls the qubit 's 
Hamiltonian [l5(. This model is relevant to a Josephson 
qubit at an optimal point, investigated in recent experi- 
ments . To emphasize specific features of this model, 
we first recall the description of the dissipative qubit dy- 
namics when the noise is either short-correlated or Gaus- 
sian longitudinal. Then we motivate model (01 and show 
that it can describe the effect of both longitudinal and 
transverse low-frequency noise. We discuss the statistics 
of X and V, and then analyze dephasing in model (J2J. 

For weak short-correlated noise the dynamics is de- 
scribed by the Bloch equations, regardless of the noise 
statistics. The weak dissipative effects from many uncor- 
related time intervals ~ r c (the correlation time) accu- 
mulate and the effect of both longitudinal and transverse 
noise can be described by the markovian Bloch equa- 
tions. The dissipative rates are given by the golden rule: 
the ovrelaxation rate 1/Ti = Su(w = e)/2 and the rate 
of dephasing (decay of a x . y ) 1/T 2 = 1/(2T 1 ) + 1/T|, the 
'pure' dephasing rate 1/T 2 * = Sy{uJ = 0)/2 being dom- 
inated by low to. This approach applies for weak noise 
with a short correlation time t c -C T±, T 2 . 

One can still rely on the lowest order of the cumulant 
expansion (but beyond the golden rule) for Gaussian lon- 
gitudinal noise V (U = 0), even for long correlations. The 
coherence (er_(t)) (here rr_ = (cr x — ia y )/2) decays then 
as 
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where Sy = ^(([V'(t), V(0)]+)) w is the noise power. 
For instance, for a linear coupling, V(X) = aX, to 
a Gaussian-distributed 1/f (flicker) noise X, when T 2 * 
defined above vanishes, one finds the dephasing law 



exp(-a 2 A f 2 i 2 |ln(w ir i)|/27r) (cf. Refs. 
the magnitude of the noise, 

S x {u)=Xf/\u\, 



Here Xf sets 



(4) 
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and may depend on the external conditions, such as tem- 
perature. The infra-red cutoff Wi r may be set, and con- 
trolled, by the details of an experiment. For instance, 
when a measurement of dephasing, performed over a 
short time i, is averaged over many runs, the fluctu- 
ations with frequencies down to the inverse of the to- 
tal signal acquisition time contribute to the phase ran- 
domization 7] (this averaging improves the accuracy and 
is needed to monitor the oscillations of expectation val- 
ues 0, IU 0> H)- In an echo- type experiment the con- 
tribution of slow fluctuations can be suppressed . The 
controlled frequency of the compensating echo pulses sets 
uiir (in Ref. Wi r ~ 1/t but it can be varied). Here we 
focus on the case w^t > 1. 

Typically, in the mentioned regime of linear coupling to 
1/ / noise dephasing is strong. To increase the coherence 
time one tunes the system to a point where the linear 
term vanishes, a = |2(. The effect of the remaining 
quadratic coupling is described by the model 10) ■ More- 
over, this model can also account for low-frequency trans- 
verse fluctuations U\{. Indeed, in the adiabatic approx- 
imation we diagonalize to — i yj(e + V) 2 + Uf f cr 2 ~ 
— |(e + Kff) c z , where V e g = V ' +Uf { /(2e). Hence as long 
as the relaxation due to the resonant (uo ~ e) part of U 
is negligible, it is sufficient to analyze the model J3J. 

This analysis requires an account of higher orders, and 
thus the knowledge not only of SV(w) but of the full 
statistics of V. The latter depends on the statistics of the 
basic quantity X as well as on the qubit-noise coupling. 
In this letter we consider a Gaussian-distributed X (see 
below), but due to the nonlinear coupling in Eq. J3J) the 
qubit is subject to non-Gaussian fluctuations V = XX 2 . 
In combination with long correlations of 1// noise, this 
shows that further analysis is needed. 

The statistics of the low-frequency fluctuations of the 
basic quantity X deserves further discussion. The as- 
sumed Gaussian statistics is generic for noise produced 
by many microscopic modes, due to the central limit the- 
orem and regardless of the noise mechanism. As for the 
flicker ('1//') noise in mesoscopic circuits, it may be dom- 
inated by only a few modes (bistable fluctuators) or con- 
tain comparable contributions of many of those, depend- 
ing on the sample, and one may expect non-Gaussian 
resp. Gaussian noise. This was demonstrated in experi- 
ments at kHz- and lower frequencies 8] . While the noise 
at higher frequencies MHz-GHz, relevant for the dephas- 
ing of qubits, is less explored, recent data suggest that it 
may have the same nature 0, 13 • Here we focus on the 
analysis of the influence of Gaussian noise X; the effect of 
a few bistable systems was discussed recently in Ref. 
Our analysis on one hand, indicates interesting features 
of decoherence for quadratic coupling; on the other hand, 
one can consider the qubit as a probe of the noise, thus 
our findings may help in identifying the noise mechanism. 

Analysis of dephasing in model To be definite, 

we assume that A > 0; the sign change has no effect 



on the dephasing laws, but reverses the phase accumu- 
lated due to the nonzero average of V. We follow the 
evolution of the off-diagonal entry of the qubit's density 
matrix (<7_(i)) = (S^cr-S), with the evolution opera- 
tor S = Texp(-%f*Vcr z dt'), which yields (<T-(t)) = 
p(i) e te*/ft( CT _(0)), where 

P(t) = (f exp f I J* Vdt'^j Texp Q J* Vdt'^j ^ , (5) 

with averaging over noise realizations. T and T denote 
time resp. reverse-time ordering; their combination in 
Eq. (J5J corresponds to the Keldysh-time ordering. Note 
the same signs in Texp and Texp, selected by the oper- 
ator <7_. 

Gaussian approximation. In the lowest-order pertur- 
bative analysis one can use Eq. © [13. For Ohmic fluc- 
tuations of X, with noise power Sx(&) — tuj coth(w/2T), 
one finds S x ^(u) = (r 2 /3tt)u(lu 2 + 4tt 2 T 2 ) coth(w/2T) 
and Eq. J2J yields for weak noise XrT <C 1 the exponen- 
tial decay with rate 
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This dephasing is stronger suppressed by cooling com- 
pared to the case of linear coupling, when 1/T 2 * ~ T. 

For 1// noise one finds S x i = (4/w)Xf In |w/ii;i r |/|w| 
and P(t) = exp(— [rfHn(wi r t)/7r] 2 ), where 

T f = XXf . (7) 

Higher orders: results and discussion. Below we ana- 
lyze effects beyond this Gaussian approximation. For the 
Ohmic noise we confirm that Eq. © holds at all relevant 
times, yielding the decay rate ©■ In contrast, for the 
1// noise corrections are strong: 



\P(t)\ = 



U)\ T t 

r ( t/2 



-1/4 



T{t « 1 , (8) 
T f t > 1 . (9) 



At short times we find a very slow decay with time scale 
(FfLr) -1 , where Lr = ln(Tf/uii T ) and we assume Lr ^> I 
(our qualitative results persist down to Lr ~ 1): At very 
short times it reduces to 1 — [Tft ln(wi r t)/7r] 2 and coin- 
cides with the result of the Gaussian approximation. This 
slow initial decay may be advantageous for application of 
quantum error correction. At longer times t ~ (TfLr)^ 1 



the decay crosses over to a power law cx 1/yt. Finally, 
at the (parametrically longer) time t ~ TjT 1 the decay 
becomes exponential ©, due to the high-frequency con- 
tribution. For large Lr > 1 a substantial decay occurs 
already in the range of Eq. JHJ . 

This unusual decay law translates into a peculiar line 
shape of the transverse spin susceptibility Xct i ( w ) = 
P(oj - e)/2 (at T < s), shown in Fig. [2 Eq. © gives 
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FIG. 1: Line shape x" x (w), in units of l/(FfLr) and as a 
function of (u) — e)/(T{Lr)- Eq. © sets a singular universal 
shape, oc (cu — e)~ 1//2 . The high-frequency contribution ||§1 |12H 
washes out this singularity on scale Tf (I/jT 1 in these units) 
and shifts the peak. The dashed line shows the result of the 
Gaussian approximation F2 combined with the phase shift F\ . 



a singular peak oc (oj — e) 1 ^ 2 , and the term 10 washes 



it out on the scale Ff, setting the peak height ~ l/yTf. 
Apart from the dephasing, we find a phase contribution 
which shifts the peak by Tf ln(o; c /rf)/7r (if this logarithm 
exceeds InLr; w c is the unltraviolet cutoff). 

Derivation. The Keldysh-time-ordered exponent © 
may be expanded into the linked-cluster series: 



P(t) = exp (f^) 



(10) 



F n representing the contribution of all connected dia- 
grams of order n in the perturbation V (the Gaussian 
approximation © neglects F n> 2). For the model J5J they 
are visualized in Fig- EJt, with a single cluster in each or- 
der. Here the solid lines represent the bare 2x2 Keldysh 
Green functions of the bath D = -iQ^Xtt)X(0)), 
whose Keldysh and retarded components [121 1 1 3| are re- 
lated to the noise power and the response function of the 
bath, respectively: D K = —2iSx, D n = ~X X - Each ver- 
tex contributes a factor A and integration over the time 
interval (0;t). Thus we find in the n-th order: 



F n (t) = 



(-A)" 



tr 



dtx / dt 



D{t l -t 2 )lD{t 2 -h)l... 
= - — ^— tr J dhj\duj2 ■ ■ ■ 

b{u}i)5 t {ui-U2)b{u2)5t{u2-W3) . . . ,(11) 

where 5t(u>) = sm(u)t/2)/(noj). The vertices contribute 
the identity matrices 1 in the Keldysh space, rather than 
the familiar f z , as a result of the sign structure in Eq. J^J ; 
such a perturbation is called 'quantum' [l^. For the cou- 
pling the averaging in JSJ reduces to a Gaussian in- 
tegral and thus to a determinant of an integral operator. 
From this standpoint our calculation may be viewed as 
determinant regularization. 

Tractable regimes. The series (jlUf) . <|11H should be eval- 
uated for each particular time t. In general, it is difficult 
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FIG. 2: a. Linked-cluster expansion. The factor 8t(Au>) = 
sin(Aw/2)/(7rAw) at each vertex violates the frequency con- 
servation, b. 'High'- and 'low'-frequency regions dominating 
the integration in the cluster diagrams. 



due to its complex structure: the propagators D are di- 
agonal in frequency and the vertices are diagonal in time. 
However, one can evaluate the dephasing in certain lim- 
its: (a) at long times (the exact condition to be speci- 
fied) St ~ <5, and frequency is conserved at the vertices; 
(b) at short times the lines become time-indepednent, 
D(ti — ti+i) m D(At — 0). These ideas are used below 
for the case of 1// noise, when D(At = 0) diverges, to 
find the short- and long-time behavior of P(t). 

1// noise: (a) higher frequencies. For the 1// noise 
it is technically more convenient to discuss the contribu- 
tions of two frequency ranges at given t (see Fig. Ep), 
rather than the short- and long-time limits. 

We begin with higher frequencies: In the integral llljl 
the adjacent frequencies may differ by ~ l/t due to the 
vertex factors S t . In the (i-dependent) range \oj\ ^> l/t 
such a shift does not change much the propagators (noise 
power), and St(Auj) ~ 5(Au>) at each vertex. Thus we 
find the contribution of high frequencies: 



\nP hi (t) = -t f — In (1 - 2iXS x 
J~i/t 2n 



H) • (12) 



The contribution of frequencies of order l/t is only es- 
timated by this expression. However, at long times 
Tft ^> 1, when Eq. I|12|) dominates P(t) (as we find 
below), this contribution is negligible, and we obtain 
Eq. ©: ln|P hf (t)| = -T^t, where 



Too — Tf 



dx 

— In 

2tt 



I -X) -I 



(13) 



As for the phase shift, to the logarithmic accuracy in the 
limit ln(r f t) > 1 we find Im lnP hf (t) = r f t ln(T f t)/ir. 

1/f noise: (b) lower frequencies. For the analysis of 
the contribution of low frequencies <C l/t one may 
replace the vertex factor St(Au) by its value, t/2%, at 
Auj = 0, and the correlations between frequencies of ad- 
jacent lines are irrelevant. Thus the scries P| gives 



F =- 



1 I 2it 



Tfduj 



— itTf In 

7T LOwt 



lnP lf (t) = -- In ( 1 - - itT { In — 

2 V 7T UJi r t 



,(14) 



(15) 
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Again, the contribution of frequencies close to the upper 
limit, uj ~ 1/t, is not given reliably by Eq. (|15|l . but it is 
negligible to the logarithmic accuracy at ln(l/Wi r t) ^S> 1. 

The low and high (and intermediate) frequencies con- 
tribute at all times, but at short times ln|P(t)| is dom- 
inated by low frequencies and at long times by high fre- 
quencies. These leading terms yield Eqs. (JSJ), JO}. (The 
susceptibility XX (w) enters the analysis, but in the final 
results we assumed a regular and hence negligible x-) 

Ohmic noise. Let us apply the developed approach 
to equilibrium thermal noise with quadratic longitudinal 
coupling, considering the example of an Ohmic bath. Its 
low- frequency spectral density x^-( w ) — i(D R — D A )/2 — 
no is related by the fluctuation-dissipation theorem to 
the noise power: D K = (D R - D A ) coth(w/2T). The 
variation scale of the latter is set by temperature, l/r c ~ 
T. Hence at times t> 1/T changing the frequency by ~ 
1 jt has little effect on D K , and one may use the long-time 
approach, with frequency conservation in the diagrams 
(analysis of D R / A does not change this conclusion). 

We consider the case of weak-coupling and begin 
with the lowest order F 2 . Evaluation of this Gaus- 
sian cluster involves tr D 2 = (D R ) 2 + (D A ) 2 + (D R - 
L> A ) 2 [coth 2 (w/2T) - 1]. The frequency integral of the 
first two terms (analytic in the upper/lower half-plane) 
vanishes, and we find exponential decay with rate Ipyl at 
i > 1/T (inch t ~ T|). 

As expected for weak noise, the higher orders provide 
only a small correction to Eq. ©, of order (\rui c ) 2 . In- 
deed, only the even orders contribute to dephasing, for 
which trD 2k = (D R ) 2k + (D A ) 2k + terms localized at 
u> < T. One finds that it is sufficient to evaluate the 
second-order contribution provided XD(lu = 0) <C 1. 
In the Ohmic case this weak-coupling condition reads 
ArT, Xru> c <C 1. 

Higher (non-Gaussian) orders as screening. The ex- 
ponential decay law appears surprising for 1// noise, 
with long-time correlations. In fact, it develops due to 
the screening of the long-time (low-w) fluctuations, in- 
teracting via the term similar to the screening of 
interaction in the Coulomb gas. Moreover, the calcula- 
tion of P{t) parallels that of the correlation energy of 
the Coulomb gas |l4J. The result can be found from 
the lowest-order cluster F 2 with the solid lines replaced 
by the properly renormalized (screened) Keldysh prop- 
agators D. Diagrammatically this appears natural: the 
diagrams F n>2 in Fig. [2i may be viewed as F 2 with ad- 
ditional vertices on the lines. 

As a result of the screening, qualitatively, the 1// di- 
vergence is cut off at u> ~ Tf, and the white low-w noise 
produces the exponential decay. We omit the detailed 
discussion but mention a specific property of the screened 
T>: unlike the bare D it is t-dependent (due to the t- 
dependent vertices). In other words, the screening sets 
in gradually, with the mentioned saturation at t > T^ 1 . 



Preparation effects. So far we worked under the as- 
sumption that at t = the bath and spin were disentan- 
gled and the bath was prepared in the thermal state of 
-ffbath- To see, if it can change our conclusions, we con- 
sidered a typical experiment, in which one monitors the 
decay of a superposition of the spin states in the pres- 
ence of permanent spin-bath coupling. We found that 
for 1/ / noise with a non-singular low- frequency suscepti- 
bility x( w ) our conculsions about the dephasing persist. 

In conclusion, we have analyzed the decay laws of co- 
herence of a qubit coupled quadratically to the environ- 
ment. We have shown that higher-order effects become 
important for certain noise spectra (notably, for 1//), 
and found the dephasing times I©-© in various regimes. 
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